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Abstract 

A fracture mechanics analysis based on strain energy 
release rates is used to study the effect of stitching in 
debonded sandwich beam configurations. Finite 
elements are used to model the configurations. The 
stitches were modeled as discrete nonlinear spring 
elements with a compliance determined by experiment. 
The constitutive models were developed using the 
results of flatwise tension tests from sandwich material 
rather than monolithic material. The analyses show that 
increasing stitch stiffness, stitch density and debond 
length decrease strain energy release rates for a fixed 
applied load. 

Introduction 

Composite sandwich construction is very attractive for 
use in aerospace structures (Figure 1). Composite 
sandwich structures are more cost effective to 
manufacture than traditional skin-stringer structures and 
have better fatigue resistance. Recent studies show that 
sandwich structure has the potential to be the lightest 
and lowest cost manufacturing concept for fuselage 
structure [1], While these configurations are efficient, 
sandwich panels can be easily damaged by foreign 
object impacts. This critical deficiency has prevented 
the use of sandwich construction in most primary 
aerospace structures. 

Debonds along the skin-core interfaces are considered 
limiting to a sandwich structure since the skin is free to 
“peel away” from the core resulting in catastrophic 
failure. This limitation is enough to discourage most 
aerospace manufacturers from using sandwich 
construction in primary structures. To overcome these 
limitations a method of mechanically fastening the skin 
to the core is needed. This method must be inexpensive 
and should not add a significant amount of weight to 
the sandwich structure. 

Stitching the skin and core is proposed here as a joining 
methodology. Stitched warp-knit textile composite 
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materials are currently being considered for use in 
primary aerospace structures and as facesheets for 
sandwich structures. Structures manufactured from 
these materials offer advantages in manufacturability and 
damage tolerance over conventional composite and 
metallic structures. Manufacturability is improved 
because large sections of dry carbon textile preform and 
core can be assembled and stitched with Kevlar threads 
near net shape before the epoxy resin is introduced/ 
Improvements in damage tolerance are observed because 
the Kevlar stitches tend to prevent propagation of 
debonds and delaminations that may be caused by in- 
plane and out-of-plane loadings. 

Figure 1 shows a representative section of composite 
sandwich structure where stitches, in addition to the 
resin interface, are used to attach the facesheets to the 
core. This structure may be subjected to a variety of 
loads as shown in the figure. Figure 2 shows a 
simplified configuration that may be used to gain 
insight into the behavior of more complex 

configurations. 

The objective of this paper is to study the effect of 
stitches on the response of a sandwich debond 
configuration. The parameters studied are the effects of 
stitch stiffness and stitching density, core stiffness and 
debond length. New constitutive models for the 
stitches have been determined from flatwise tensile tests 
of a sandwich material and are implemented in the 
present analyses. A debond between the top face sheet 
and the core is assumed and strain energy release rates at 
the debond front are used to evaluate the effects of 
several parameters using finite element analyses and the 
virtual crack closure technique (VCCT). 

Debonded Sandwich Beam Configuration 

A debonded sandwich beam configuration subjected to 
tensile loading is shown in Figure 2(b). The stitched 
sandwich beam structure is assumed to have Kevlar 
stitches, foam core, and carbon facesheets as shown in 
the figure. The parameters of this configuration are the 
length of the specimen, L, the thickness of the 
facesheets, //, the thickness of the core, t c , and the 
spacing of the stitches, s. The local details of a debond 
of length, a , propagating from one end of the 
configuration including the loading of the stitches are 
shown in the figure. 
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Facesheet material and thicknesses that are representative 
of a hybrid IM7/3501-6 and AS4/3501-6 warp-knit 
fabric are considered. 3 The material consists of IM7 
yams in the axial direction and AS4 yams in the off- 
axis directions. Each stack of material is assumed to be 
oriented with its primary axis in the ^-direction and 
having a thickness of 0.055 in. The overall density of 
the fabric expressed in areal weight is 0.577 oz./ft. 2 , 
1.21 oz./ft/ and 0.651 oz./ft/ for the forty-five, zero 
and ninety degree plies, respectively. The equivalent 
laminate stacking sequence of each stack of material is 
(45/-45/0/90/0/-45/45) ns where n=2 for both of the 
facesheets. 

Finite Element Analysis 

Figure 1 shows a stitched sandwich shell-type structure. 
Three-dimensional modeling and analysis of this 
complex configuration requires a large finite element 
model. However, considerable insight into the behavior 
of such a complicated configuration can be obtained by 
studying simpler configurations such as the debonded 
sandwich beam configuration shown in Figure 2 while 
reducing modeling complexity. The debonded 
sandwich beam configuration was modeled using the 
ABAQUS finite element code. 4 

The method of analysis implemented here involves the 
use of plane strain continuum elements to model the 
configuration and nonlinear spring elements with 
experimentally determined compliances to model the 
stitches. This modeling technique allows an 
experimentally determined load vs. deflection behavior 
to be considered for each stitch. The stitch response 
includes the local effects of the stitch debonding from 
the laminate in addition to nonlinearity of the laminate, 
core and stitch materials. 5 

The specimen configurational parameters including 
facesheet and core thicknesses and stitch spacing are 
shown in Figure 2(a) and Table 1. Debond lengths, a, 
between 0.250 in. (a/if/* 2.27) and 2.00 in. (a/tj= 18.2) 
are considered for the unstitched and stitched 
configurations with a length, 1=7.00 in. and an applied 
load, />= 100 lb. 

Material Properties 

In these analyses, the laminates are assumed to be 
homogeneous with axial properties determined 
experimentally and all others estimated using the 
equivalent stacking sequence and classical lamination 
theory as 


Eh=T1.7 Msi p.12-2.50 Msi \)i 2 -0.40 
E22-5.14 Msi |Xi3— 1 .77 Msi \>i 3 — 0.30 
E 33 -I .79 Msi JJ.23~0.88 Msi \) 2 3=0.30 

where (ij— 1,2,3) are the Young’s moduli, 

shear moduli, and Poisson’s ratio, respectively, and the 
subscripts 1,2,3 represent the fiber, transverse and out- 
of-plane directions, respectively. 

Three cases of core properties are considered. The first 
two cases consider a foam core approximated as an 
isotropic continuum with elastic properties (£ c , v) and 
densities (p) taken from vendor data as 

E c =5220 psi u^O.385 p l =2.0 Ib./ft* 

E c 2 = 13340 psi \} 2 =0.385 p 2 =4.68 lb. /ft 3 

The third case considers a core with the same 
orthotropic properties as the facesheets (E c 3 =E /ace , hech ). 
Core moduli and stitch spacings for all of the debonded 
sandwich beam configurations considered in this paper 
are shown in Table 1 . The three cases represent soft and 
medium density foam cores and the extreme case of a 
core having the same properties as the carbon/epoxy 
facesheets. 

Strain Energy Release Rates 

The configuration was modeled with the ABAQUS 
finite element code using 8-noded quadratic two- 
dimensional (plane strain) elements. 4 A representation of 
the 8-noded plane strain elements near a debond front 
with rectangular grid type modeling is shown in Figure 
3 with elements having length of 0.0625 in. The 
elements are assumed to have the same length, A, ahead 
of and behind the debond front (as shown in Figure 
3(b)). The G-values can be calculated using the nodal 
forces (F x> F : ) and displacements (u, w) near the debond 
front and the increment of new debond area created 
using the VCCT as (see Figure 3) 6 ' 7 

Mode-1 components: 

Gi = ~J^[ F zi( w P ~ w p‘) + F zj( w t ~ )] ^ 

Mode-11 components: 

G " = ~2l[ Fjt '(“ p ~“ p ') + F *j 

with 

G Totai = Gy + G n . (3) 
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Configuration 

Facesheet 
Thickness 
tr (in.) 

Core 

Thickness 
U (in ) 

Stitch 
Spacing 
s (in.) 

First Stitch 
Location 
(in.) 

Core 
Modulus 
E c (psi) 

N-i 

0.110 

0.500 

N/A 

N/A 

5220 

N-2 

0.110 

0.500 

N/A 

N/A 

13340 

N-3 

0.110 

0.500 

N/A j 

N/A 

E„ = En 

S4-I 

0.110 

0.500 

0.250 

0.00 

5220 

S4-2 

0 110 

0.500 

0.250 

0.00 

13340 

S4-3 

0.110 

0.500 

0.250 

0.00 

fcq Mj 

S2-1 

0.110 

0.500 

0.500 

0.00 

5220 

S2-2 

0.110 

0.500 

0.500 

0.00 

13340 

S2-3 

0.110 

0.500 

0.500 

0.00 

p core p laccsficcii 

— £-ti 

S41-1 

0.110 

0.500 

0.250 

0.250 

5220 

5TT3 1 

5TR5 

oro 

515(5 

0.250 

13340 

S41-3 

0.110 

0.500 

0.250 

0.250 

corc_ r- laccsflects 
Ml ~ Mi 

S21-1 

0.110 

0.500 

0.500 

0.500 

5220 

S21-2 

oTTo 

i 0.500 

0.500 

0.500 

13340 

S21-3 

0.110 

0.500 

0.500 

0.500 

p COft p lilCCShCCU 

Mi Mi 


The strain energy release rates of the debonded sandwich 
beam configurations are calculated using equations (1)- 
(3) with the nodal forces and displacements transformed 
into a coordinate system that is normal to the debond 
front in a geometrically nonlinear analysis. For cases 
with a stitch present at the debond front nodal location, 
the stitch force is subtracted from the computed force at 
the debond front in the VCCT calculations consistent 
with the calculations in reference 7 for pressure loaded 
debond faces. 

This configuration may exhibit an oscillatory stress 
field very near the debond front and demonstrates 
nonconvergence of the individual modes of the strain 
energy release rate. The computed values of G { and Gu 
do not converge with decreasing element size, however, 
the sum of the modes does converge and is given as 
G T o,«i in Eq. (3). 78 

Modeling Stitches 

The configurations have been analyzed with a 
geometrically nonlinear finite element analysis within 
the ABAQUS finite element code. Using similar 
techniques discussed in references (9-1 1), the stitches are 
modeled as nonlinear spring elements rather than as spar 
or beam elements. 

The spring elements are assumed to have only axial 
stiffness, K m mi . Unlike the analyses discussed in 

references (9-11), spring elements ahead of and behind 
the debond front are considered because of the through- 
the-thickness deformation of the compliant foam core. 
The nonlinear spring elements are evenly spaced along 
the debonded length of the model. The stitch location 
spacing presented in Table 1 coincides with element 
comer nodal locations. 


Compliance curves for both axial and shear behavior of 
the stitches embedded within a laminated carbon/epoxy 
textile composite with material properties similar to 
those of the sandwich facesheets were developed in 
reference 5 using flatwise tension and double lap shear 
tests. Similar curves have been developed for a Kevlar 
stitch embedded within a sandwich material. 
Constitutive behavior of the stitch was determined from 
flatwise tensile tests of a sandwich composite with a 
Teflon insert separating the facesheets from the core, 
thus allowing the nonlinear response of the stitches to 
be determined. 

Two flatwise tensile specimens were tested and curves 
representing their load-deflection behavior (stiffness) are 
labeled as Specimen I and Specimen 2 in Figure 4. A 
curve representing the average load-deflection behavior 
of the two specimens along with curves representing 
hypothetical upper and lower bounds are also shown in 
the figure. The upper and lower bounds on likely stitch 
stiffnesses are assumed to be 150% (stiff stitch values) 
and 50% (compliant stitch values) of the average 
stiffness, respectively. These compliance curves 
represent the net behavior of the stitch due to stitch, 
core and facesheet material nonlinearity and stitch 
debonding. 5 

A piecewise linear representation of this data is used in 
the finite element model. Examination of the flatwise 
tensile test data revealed that the axial stiffness of the 
stitches increases from approximately 3300 lb. /in. 
initially to 5500 lb./in. near failure. Note that failure of 
the stitches occurs at an average load of 60 lb. per 
stitch. These stiffnesses and failure loads are used for 
the characterizations in this paper. 

Each stitch is modeled with a single nonlinear spring 
element joining the inner surfaces of the upper and 
lower facesheets. Thus, compatibility between the 
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stitches and core is not maintained, nor is contact 
between the stitches and the core considered. The effect 
of this simplification may be minor when considering 
compliant cores, but should become more important as 
the stiffness of the core increases. 

Modeling the Facesheet to Core Contact Problem 
Closure of the debonded faces may occur once the 
debond is of sufficient length because of the stitch 
forces. In the finite element analysis, contact of the 
faces is allowed, while interpenetration of the faces is 
not. In references (9-11), interpenetration of the faces 
was prevented by adding multipoint constraints along a 
known region of interpenetration to impose the 
requirement of identical z-direction (w) displacements 
among elements in contact. 

Because of the compliance of the core in the sandwich 
structure, the stitches are assumed to carry load in 
compression as well as tension. In the present analysis, 
the stiffness of the stitches in compression is assumed 
to be equal to the initial stiffness of the stitches in 
tension. Given this assumption, the computed values 
of facesheet to core interpenetration typically were less 
than 0.2% of the facesheet thickness. Thus, multipoint 
constraints are not used in the present analysis. 

Results and Discussion 

Four variables are considered in this analysis and 
include stitch stiffness, stitch configuration, core 
stiffness and debond length. Three stitch stiffnesses are 
considered: the average value ffom testing and 150% 
and 50% of the average value. Five stitch 
configurations are considered: (a) unstitched, (b) four 
stitches per inch, (c) two stitches per inch, (d) four 
stitches per inch with no stitch under the load, and (e) 
two stitches per inch with no stitch under the load. 
Configurations (b) and (c) have their first stitch directly 
beneath the load while configurations (d) and (e) have 
their first stitch offset ffom the load by a distance equal 
to the stitch spacing. Three core stiffnesses are 
considered: a compliant core with modulus, £'=5220 
psi; an intermediate stiffness core with modulus, 
£‘=13340 psi; and a stiff core with modulus, 

£ -Ejdc'shni,. Finally, debond lengths, a, between 
0.250 in. (a/t/= 2.27) and 2.00 in. (a/tj= 18.2) are 
considered with an applied load, P, of 1 00 lb. 

As shown in Figure 5, the unstitched configurations 
exhibit increasing G To ui with increasing debond length, 
a, for all three moduli (configurations N-I, N-2 and N-3 
in Table 1). However, increasing core modulus 
decreases Great even for large debond lengths. 


In the discussion of the stitched configurations shown 
in Figures 6 through 10, results are given for (a) stiff 
stitch values, (b) average stitch values and (c) compliant 
stitch values as shown in Figure 4. 

The effect of stitching is shown in Figures 6(a) through 
6(c) for configurations with four stitches per inch ( S4-I , 
S4-2 and S4-3). The slope of the G To ui curve for 
stitched configurations with foam core ( S4-1 and S4-2 ) 
decreases for all stitch stiffnesses and all debond lengths 
because of the load transfer of the stitches. The slope of 
the Great curve does not begin to decrease for the 
configurations with the stiffest core (S4-3) until the 
third stitch begins to carry load. The G’i 0 , a: continues to 
decrease with both the number of loaded stitches and 
the distance from the debond front to each of the load- 
bearing stitches. For long debonds (atf>9 ()9). the 
strain energy release rate approaches near-zero values in 
the stitched configurations. Note that the G-values are 
functions of stitch stiffness but even the compliant 
stitches have a profound effect on the G-values as the 
debond length increases. 

The effect of a lower stitching density (2 stitches per 
inch) is shown in Figures 7(a) through 7(c). As with 
the S4 configurations, the S2 configurations (S2-L S2-2 
and S2-3) show a general decrease in Great with 
increasing a/tj. The reduction in strain energv release 
rate requires a longer length of debond growth because 
there are fewer load bearing stitches per inch of debond 
length. G-values are somewhat larger for the 
configurations with lower stitching density than for the 
configurations with four stitches per inch but this effect 
is not as significant as that of the stiffness of the 
individual stitches themselves (e g. compare Figure 7(a) 
and 7(c)). 

The effect of stitch placement is considered next. Note 
that a load bearing stitch is present at x=0 in. 
(immediately below the point of load application) in 
both the S4 and S2 configurations. Figures 8(a) through 
8(c) show Great plotted against debond length for 
configurations with four stitches per inch and missing 
the stitch immediately beneath the load (547-/, S4I-2 
and 547-3). Compared with configurations 54 (Figure 
6), the maximum values of Great for configurations 54 7 
(Figure 8) have increased significantly. Although the 
effect of the stitches is delayed and somewhat reduced, 
the analyses for long debond lengths {a/tj> 9.09) shows 
that Great decreases toward zero as the debond grows 
even with sparse stitching. 

Finally, the extreme case of two stitches per inch with 
no stitch under the load (527-7, 527-2 and 527-3) is 


4 

American Institute of Aeronautics and Astronautics 



considered in Figures 9(a) through 9(c). Note that the 
maximum value on the vertical scale has increased from 
2.0 in. -lb. /in. 2 to 5.0 in.-lb./in. 2 Again, the stitches are 
shown to reduce the G-values to near zero for long 
debond lengths, (a/tj= 18.2). 

The decrease in G Total for longer debond lengths is the 
result of the stitches resisting debond opening (reducing 
the opening force in the region near the debond front.) 
This effect increases with increasing distance between 
the debond front and each load bearing stitch and also 
with increasing numbers of load bearing stitches. 

Concluding Remarks 

The effect of stitches on the total strain energy release 
rate for a debonded sandwich beam configuration with 
an initial debond was studied. Plane strain finite 
elements were used to model the configuration and the 
virtual crack closure technique was used to calculate the 
strain energy release rates. Using this fracture 
mechanics approach, a debond between the skin and 
stiffener flange was assumed. The debond growth 
between the facesheet and core was assumed to be self- 
similar and continuous along the length of the interface. 
The stitches were modeled as discrete nonlinear spring 
elements with their compliance determined by 
experiment. Interaction between the stitches and core 
was not considered. 

For the unstitched configuration, Gw decreases with 
increasing core modulus. The effect of stitching was 
most significant for compliant configurations. As the 
debond length increases, additional stitches begin to 
carry load resulting in a decrease in Gtoui • For long 
debonds, the stitches may produce enough compressive 
force at the debond front to close the debond and reduce 
G Total to near-zero values. The stitching density, stitch 
stiffness and proximity of the first stitch to the point of 
load application effect the rate at which Gt 0 ui decreases. 
Configurations with a dense stitch pattern and stiff 
stitches decrease Gtotai more rapidly than configurations 
with sparse or compliant stitching. 
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Figure 1 . Stitched sandwich concept. 




(b) Deformed configuration and boundary conditions 



(c) Typical finite element mesh 


Figure 2. Debonded sandwich beam configuration. 
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G Total 

(in.-lb./in. 2 ) 



Figure 5. Strain energy release rate for unstitched configuration. 




(a) Stiff stitch values 


fb) Average stitch values 


Glotai 

(in.-lbVin. 2 ) 



Debond Length, a/t, 

(c) Compliant stitch values 


Figure 6. Strain energy release rate for configuration with 4 stitches per inch. 


8 

American Institute of Aeronautics and Astronautics 





0.00 5.00 10 00 15.00 20.00 


Debond Length, a/t, 

(a) Stiff stitch values 



000 500 10.00 1500 2000 


Debond Length, a/t, 


(b) Average stitch values 


G 


Total 


(in.-lb./in. 2 ) 



Debond Length, a/t, 


(c) Compliant stitch values 

Figure 7. Strain energy release rate for configuration with 2 stitches per inch. 
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Debond Length, a/t, 

(a) Stiff stitch values 


GlotaJ 

(in.-lb./in. 



Debond Length, a/t, 

(b) Average stitch values 



Debond Length, a/t, 


(c) Compliant stitch values 

Figure 8. Strain energy release rate for configuration with 4 stitches per inch (no stitch under the loading 
point). 
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G Total 

(in.-lb./in. 2 ) 


500 


4 00 


300 


2 00 


1 00 
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10 00 
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Debond Length, a/t, 

(a) Stiff stitch values 



0.00 5.00 10 00 1 5.00 20.00 

Debond Length, a/% 

(b) Average stitch values 



0 00 5 00 10 00 1 5.00 20 00 

Debond Length, a/t r 

(c ) Compliant stitch values 

Figure 9. Strain energy release rate for configuration with 2 stitches per inch (no stitch under the loading 
point). 
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